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ABSTRACT 


This thesis develops a methodology for evaluating mine actuation data. 
Computer models are developed to analyze actuation data obtained from the Mine 
Warfare Command by fitting various types of actuation curves to the data. For 
each actuation curve type, maximum likelihood estimates are used to determine 
those parameters resulting in the greatest probability of obtaining the observed 


data. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not have 
been exercised for all cases of interest. While every effort has been made, within the time available, 
to ensure that the programs are free of computational and logic errors, they cannot be considered 
validated. Any application of these programs without additional verification is at the risk of the 


user. 
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EXECUTIVE SUMMARY 


In mine warfare, a commander’s ability to make the correct decisions while conducting mine 
countermeasure (MCM) operations could mean the difference between life and death for those going 
in harm’s way. Over the years, a number of computer models have been developed to assist 
commanders in making those difficult decisions. However, the outputs generated by these useful 
tactical decision aids are only as good as the tactical parameters they used for inputs. The analytical 
techniques used for estimating the input parameters are crucial to the success of MCM operations. 
Two of the more commonly required input parameters are actuation width (4) and actuation 
probability (B). 

This thesis develops a methodology for the evaluation of mine actuation data. The analytical 
approach used produces maximum likelihood estimates of A and B. Separate analysis methods are 
also developed to test the symmetry/asymmetry of actuation data and to estimate the mean and 
standard deviation of a ship’s navigational error. 

The resulting methodology is used to analyze mine actuation data provided by the Mine 
Warfare Command The data were obtained during an MCM exercise in a realistic scenario. The 
analysis concludes that this data set is not statistically inconsistent with the assumption that the 
observed data was generated from a symmetric actuation curve. The mean and standard deviation 
of the minesweeper’s navigational error are estimated to be 19.74 yards and 159.2 yards 
respectively. The data was fit to three different types of generalized actuation curves with the 


following results: 


X11 


Rectangular Symmetrical | Rectangular Asymmetrical {| Washburn Actuation 
Actuation Curve Actuation Curve Curve 
2250 2125 2250 


3388 3333 3400 
295x.10" 275 x10" 2.25x 10 





Note: C in the table above is a parameter for the Washburn Actuation Curve. 
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I. INTRODUCTION 
A. BRIEF HISTORY 
The underwater mine, it can be said, came of age with the 


First World War, it matured during the Second World War, and now its 
future as one the principal weapons of attrition and defence seems assured. 


[Ref. 1] 

Throughout history naval mines have been effectively used to deny the enemy 
presence within coastal waters and to ensure friendly control of strategic waterways. The 
use of naval mines as a weapon is a concept which dates back to 668 A. D. In more recent 
times, naval mines have played a significant role in every major war involving the United 
States. 

During World War I, the Germans successfully employed defensive and offensive 
mine warfare against the British. Most Aobie were the results obtained from the German 
minefields laid west of Heligoland. In the Battle of Jutland the British decided not to follow 
the German Fleet into the Heligoland Bight because of their concern over the German 
minefields. In a similar fashion, the British also benefited from the use of naval mines 
during the war. They were very successful in denying German submarines the use of the 
tactically significant Dover Straits. British minefields consisting of over 5,000 mines forced 
the Germans to discontinue the use of these straits and reach the Atlantic through the longer 
northern route, thus alleviating British merchant shipping losses.[Ref. 2] 

World War II brought about significant advances in mine technology. New | 
developments such as more sophisticated firing mechanisms, safety and delay mechanisms, 


and ship counters made mine warfare an even more serious threat. A total of over 300,000 








mines were laid by the United States and Great Britain during the war. The total loss of 
enemy shipping attributed to mining projects is estimated to be close to 2,700 ships [Ref. 3]. 
Post-war information received from senior Japanese naval officers indicated that the mining 
of Japan and the numerous losses of merchant vessels due to mines were factors which 
greatly influenced the outcome of the war. 

The delay of D-Day at Wonsan during the Korean War clearly highlights the 
significance of possessing a well-trained, well-equipped mine countermeasures (MCM) force 
capable of conducting swift and effective mine clearance operations. Wonsan served to 
unveil the U.S. Navy’s limited capabilities in shallow-water mine clearance spemtions 
During the minesweeping operations at Wonsan, the United States suffered the loss of four 
minesweepers, one fleet tug, and had five destroyers severely damaged by mine blasts. A 
massive minefield containing an estimated 3,000 mines paralyzed a fleet of 250 ships 
carrying 50,000 Marines off the coast of Wonsan for one week [Ref. 4]. After the events at 
Wonsan, Admiral Forrest P. Sherman, Chief of Naval Operations stated: 

When you can’t go where you want, when you want to, you 
haven’t got command of the sea. And command of the sea is a rock bottom 
foundation of all our war plans. We’ve been plenty submarine-conscious and 
air-conscious. Now we’re going to start getting mine-conscious--beginning 
last week.[Ref. 5] : 

The vast majority of MCM operations during the Vietnam War took place in the 
rivers of Vietnam. The North Vietnamese relied heavily on the mining of key inland 
waterways to control their use for logistical purposes. The shallow water MCM operations 


in Vietnam were aggravated by the Vietnamese use of combined gunfire and rocket attacks 














against MCM forces, making sweeping operations more difficult and dangerous than they 
were to begin with. One of the major mining campaigns during the war was the U.S. mining 
of Haipong and other North Vietnamese harbors. The results of the almost 8,000 mines laid 
were immediate and effectively stopped all ship traffic in or out of the harbors for ten 
months. 

History repeated itself during the Gulf War against Iraq. In preparation for an 
amphibious assault, coalition forces were ordered to sweep a channel in order to provide safe 
passage for a battleship to a Fire Support Area off the coast of Kuwait. During the clearance 
operation the flagship of all MCM forces in the Gulf and its anti-air warfare coverage ship 
were struck by mines. The Chief of Naval Operations, Admiral Frank Kelso best 
summarized the events when he said: 

We recently relearned some hard lessons-- how mines can 
frustrate even the most powerful of naval forces. During Operation Desert 


Storm, Iraq’s extensive minefields all but stymied a planned amphibious 
strike to liberate Kuwait.[Ref. 6] 


B. OBJECTIVE 

The main objective of this thesis is to explore the use of a fresh analytical approach 
for the evaluation of exercise mine actuation data. An analysis is performed of actuation 
data obtained during an MCM exercise in a realistic scenario. Maximum likelihood 
estimates (MLE) are used to determine those values of actuation width (A), and probability 
of actuation (B) which best represent the exercise data. A separate analysis is used to 
determine the standard deviation (go) and mean () of the navigational error which were 


needed for the MLE analysis. 











C. OUTLINE 

Chapter II introduces general concepts in mine warfare. Chapter III reviews 
preliminary analysis techniques used by the Mine Warfare Command, discusses the analysis 
methods used in this thesis, and presents the results of the analysis conducted. Chapter IV 


concludes the thesis work and comments on areas of possible future work. 











Il. BACKGROUND 


The effectiveness of the (submarine) mine has not decreased 

with the coming of the space age. So long as cargo ships cross the sea, this 

unspectacular weapon will remain a major factor in control of the approaches 

to harbors, and the shallow straits between seas.[Ref. 7] 
A. MINE WARFARE 

Due to its relative low cost and high effectiveness, the use of mines is a tactic which 
we are likely to see in future wars. Today’s global political, economic, and military 
conditions have forced the U.S. Navy to shift its attention from blue water operations to the 
littorals. And it is in the world’s littorals where the capability of conducting effective and 
thorough mine warfare becomes critical. Admiral Kelso’s viewpoint : “Effective offensive 
and defensive mine warfare underlies--literally--the success of littoral military operations 


[Ref. 8].”, is representative of the importance afforded to mine warfare by the Navy’s senior 


leadership. Mine warfare serves two main objectives: 


@ To damage or destroy enemy shipping. 
@ To deny the enemy use of certain waters, or at least hinder his operations in these 
waters by the threat presented by a minefield. 
1. Types of Mines 
Naval mines are commonly grouped into two main categories: 
a. Controlled mines 


Controlled mines are those whose detonation is electrically controlled from 








a shore site. These mines can be set to off to allow safe passage of friendly shipping, or they 
can be set to on to be used against enemy shipping: 

b. Independent mines 

Independent mines act independently once laid and do not require user 
intervention for detonation. They are equally lethal to all shipping regardless wether they 
are enemy, neutral, or friendly ships. Independent mines can be further classified by the 
position they occupy in the water. 

(1) Moored mines. This type of mine is buoyant and floats at a 
predetermined depth below the surface. It is held in position by a cable attached to an 
anchor. 

(2) Ground mines. Also known as bottom mines, these mines lay on 
the bottom of the ocean. They are very effective against shallow water shipping and pose 
a significant threat to submarines in deep water. 

(3) Drifting mines. These are mines which move freely at or near the 
surface of the ocean. They are buoyant or neutrally buoyant, and can be attached to a float 
line at a set depth beneath the surface. 

(4) Creeping mines. A type of drifting mine held below the surface 
by means of a length of wire or chain which drags along the bottom. 

(5) Oscillating mines. These mines are free floating with a 
predetermined range of upper and lower depth below the surface. They use compressed air 


or gas to maintain their position in the water. 








Furthermore, independent mines can also be categorized by the type of 
actuation mechanism which they use to detonate themselves: 

(6) Contact mines. Mines which detonate only through physical 
contact with a ship. 

(7) Magnetic mines. These mines are actuated when they detect a 
disturbance in the earth’s magnetic field, such as the one caused by a steel hull vessel. 

(8) Acoustic mines. These mines possess hydrophones which are 
tuned to detect mechanical noise made by shipboard machinery. They are actuated when 
noise within a specified frequency band is detected. 

(9) Pressure mines. These mines are actuated by sensing the pressure 
variations in the water directly underneath a moving ship. 

(10) Combination mines. These mines utilize at least two of the three 
basic influence mechanisms (magnetic, acoustic, or pressure). They are actuated when all 
its influence actuation criterion are met. 

Zi Minefields 
There are two basic kinds of minefields: 
a. Defensive 
The early minefields were basically defensive in nature. These fields were 
designed with the main objective of keeping enemy ships out and protecting friendly 
shipping. This objective is still valid for present time minefield planning. Defensive 
minefields are mainly used in friendly harbors and strategic waterways under friendly 


control. The idea is to allow safe passage of friendly vessels while denying the enemy the 











opportunity to position itself to conduct attacks, invasions, or disruptions of friendly 
operations. The mining of Wonsan during the Korean War is an excellent example of how 
effective defensive minefields can be against the enemy. 

b. Offensive 

Offensive mining is a concept which was spurred by the ability to lay mines 
from aircraft and submarines. In the early days of mine warfare any offensive mining had 
to be done from surface ships, which have a more limited ability to operate in opposed 
waters. These minefields are used to attack enemy shipping and to deny the enemy use of 
its harbors and waterways. As such, offensive minefields are laid in waters controlled by 
the enemy. During the Gulf War, Iraq forced American commanders to reconsider a planned 
amphibious assault when two U.S. warships struck mines laid by Iraq in the waters off the 
Kuwaiti coast. 

3. Mine Countermeasures 


There are two basic methods employed during the conduct of mine countermeasures: 


@ Use of special equipment to reduce the ship’s signature. 


@ Physical removal or disarmament of the mine. 


a. Signature reduction 

Modern warships have a limited capability to actively protect themselves 
against various types of mines. In the case of magnetic mines, today’s ships attempt to 
counter the threat by means of degaussing coils. Degaussing coils are intended to counteract 


the disturbances a ship constructed out of steel produces in the earth’s magnetic field. 























Against acoustic mines, engineering technological advances such as sound mounts and 
dampening materials have made it possible to build engines and machinery. which operate 
at much lower noise levels than in the past. The only countermeasure in use against pressure 
mines is simply the reduction of ship’s speed in an attempt to minimize water disturbances 
caused by a ship’s motion through the water. Moored mines are in general easier to defeat 
as they can be detected by sonar, lookouts, or helicopters. 

b. Physical removal or disarmament 

As occurred during the Gulf War, naval vessels a sometimes ordered to 
operate in waters which are known or suspected to be mined. if such cases, essential routes 
through the water are first sanitized by conducting minesweeping and minehunting 
operations. 

(1) Minesweeping. Minesweeping operations are tailored to the 
specific type of mine being swept. For influence mines, either magnetic or acoustic devices 
are used to simulate a ship’s influence field. The acoustic devices are underwater 
mechanisms which produce noises similar to shipboard engines and machinery. In principle, 
the object is to fool the mine’s acoustic sensors causing it to actuate and explode. The 
magnetic influence devices are essentially an electrically charged cable being towed from 
a ship or a helicopter. With this type of device the goal is to create an electromagnetic field 
strong enough to disturb the vertical component of the earth’s magnetic field and thus 
actuate the mine. For moored mines, a paravane device is used to support a cable at its outer 
end while holding it out at an angle to the sweeper. As the cable is towed through the water, 


cutting blades attached to it cut the mooring lines of mines laying in its path. The mines then 





surface and are destroyed by Explosive Ordnance Disposal personnel. 

(2) Minehunting. Minehunting refers to the use of sophisticated sonar 
equipment to identify minefields and accurately classify the mines within them. This is a 
tedious and extremely difficult operation whose success is not only dependent on the 
sensitivity of the sonar system, but is very much influenced by the environmental and ocean 
floor conditions. Minehunting had its origins during World War II, where one of the 
preferred methods employed to locate mines was to have mine watchers guard key 
waterways to spot the position where mines hit the water after having been laid by aircraft. 
Other possible techniques for minehunting included the use of radar to fix splash positions, 
and the use of a hydrophone network to determine splash locations by computing time 
differences between sound arrivals at various hydrophones. Since then, various types of 
minehunting sonars have been developed to find mines and neutralize them after they have 
been deployed and have reach their resting location in the water. The U.S. Navy currently 
has a class of ships designated as coastal minehunters (MHC) which are equipped with the 
latest sonar and mine neutralization systems. Research and development efforts are currently 
ongoing for a new class of autonomous underwater vehicle to be used in a minehunting role 


[Ref. 9]. 
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Ill. ANALYSIS METHODOLOGY AND RESULTS 


A. EXERCISE DESCRIPTION 

The data analyzed in this thesis were obtained during a minesweeping exercise 
conducted in 1984. Even though these data are quite old, they were used because sweeping 
procedures and mine actuation mechanisms have not significantly changed, and the data 
show a puzzling asymmetry which initially suggested that the minesweeper was considerably 
more effective on the port side than on the starboard side. This hypothesis of asymmetry 
was eventually shown not to be statistically supportable by the exercise data. 

The first step in evaluating the mine actuation data was to gain an understanding of 
the physical setting of the exercise, the information being represented by the data, and the 
data collection procedures. Appendix A contains the raw data as it was received from the 
Mine Warfare Command. 

Five influence exercise mines were laid in a simulated minefield measuring 1200 
yards in width by 3700 yards in length. The mines were deployed via small boat by an 
Explosive Ordnance Disposal team using precise navigation equipment (Mini-Ranger 
station). The navigational equipment provided an accuracy of 2-3 meters in the placement 
of the mines which allowed us to assume that any errors associated with the mine locations 
were negligible. As shown in Figure 1, the mines were laid at equal intervals along a 
diagonal line at a water depth of 15-19 fathoms. After the mines were in place, a 
minesweeper made a total of ten north/south parallel sweeps across the width of the 


minefield. For each individual mine encounter with the minesweeper, post-exercise data 


1] 








collection included actuation result (actuation or non-actuation), mine position relative to 
ship’s track (left or right), lateral range between the ship and the mine (1.e., the closest point 
of approach), and time of actuation. The actuation result and time of actuation information 
were obtained from data recording systems housed within each mine. The relative position 
and lateral range information were reconstructed after the exercise using the ship’s 


navigational charts. 











-800 -600 -400 -200 200 400 600 800 





Figure 1. Simulated Minefield. 
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B. ACTUATION CURVES 

1. Idealized Actuation Curves 

An actuation curve is a plot of the probability that a ship passing a mine at some 
lateral range (d) will cause the mine to actuate. The lateral range d is the distance between 
the mine and the ship at the closest point of approach. The ship’s track is assumed to be 
infinite in both directions and to contain no course changes. It is important to understand 
that an actuation curve is not a probability dicnsity function (PDF) nor a cumulative 
distribution function (CDF). Actuation curves simply represent the cumulative probability 
of actuation sometime along the ship’s track and under a particular set of environmental 
conditions. 

For this analysis, the exercise data was fit to three types of idealized actuation curves: 
a symmetrical rectangular actuation curve, an asymmetrical rectangular actuation curve, and 
Washburn’s actuation curve [Ref. 10]. Figure 2 shows examples of each type of actuation 


curve used in the analysis. 


C. PRELIMINARY ANALYSIS 

During the 1950's R. K. Reber conducted valuable research work on mine warfare 
issues. Much of his work and conclusions have become the basis for current analysis 
methodologies in the mine warfare field. One such methodology is the NATO Standard 
Naval Agreement (STANAG) 1142. This NATO Agreement delineates a procedure for 


estimating actuation width (A) and actuation probability (B) from mine actuation data. 
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The Mine Warfare Command uses this methodology to conduct preliminary analysis of field 


data. The analysis conducted in this thesis is aimed at exploring alternative data reduction 














procedures. 
: | 
) aid eee 
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Figure 2. Types of Actuation Curves. 
1. Empirical Actuation Curve 


The preliminary analysis performed on the exercise data (at the Mine Warfare 
Command in 1985) consisted of estimating A and B using empirical techniques from 
STANAG 1142. The underlying idea is to construct a histogram using field data to provide 


an approximation of the actual actuation curve. The histogram is generated by partitioning 
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the lateral range of the data into intervals of equal length, and determining the empirical 
probability of actuation for each interval. The probability of actuation for any given lateral 
range interval is calculated by dividing the number of mine actuations in a particular lateral 
range interval (x) by the total number of mine encounters in that particular interval (”). In 
this context, a mine encounter is a single pass between a ship and a mine. Mine encounters 
are completely characterized by the lateral range between the ship’s track and the mine, and 
the actuation outcome. Negative lateral ranges result when the mine passes on the port side 
of the ship. Each time a mine encounter takes place, the outcome is either an actuation or 
a non-actuation. 

Figures 3a and 3b show two possible empirical actuation curves based on the exercise 
data used in this thesis. In Figure 3a, a left-right symmetry is assumed by considering only 
the absolute value of the lateral range. Figure 3b shows the empirical actuation curve where 
a distinction is made between mines passing on the left or the right of the minesweeper. In 
the asymmetric curve it can be seen that the majority of the mine actuations occurred on the 
left side of the minesweeper. Upon careful review of the exercise procedures, it was 
discovered that the asymmetry in this curve was caused in large part by the design of the 
exercise. Assuming perfect navigation, the minesweeper would have seven opportunities 
to pass mines off its starboard side, 37 opportunities on the port side, and six opportunities 
with a lateral range of zero. From the reconstructed tracks, 11 mines actually passed off the 
starboard side, 37 off the port side, and 2 down centerline. This helps to explain why the 
majority of mine encounters and mine actuations were reported as having occurred on the 


left side of the minesweeper. 
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Figure 3a. Empirical Actuation Curve of Exercise 


Data Assuming Symmetry. 
Figure 3b. Empirical Actuation Curve of Exercise 


Data Assuming Non-Symmetry 








The undersampling that occurred in the starboard range bins causes some 
complications in data analysis. The actuation data in the asymmetric curve in Figure 3 
suggest that the minesweeper had very little effectiveness on its starboard side, but this 
might be explainable by the small sample sizes obtained on the starboard side. Intuition 
suggests that the empirical actuation curve should be approximately symmetric. In the 
remainder of this subsection, we develop confidence regions for the actuation proportions 
observed in each range bin. The intent of this development is to show that the assumption 
of a symmetrical actuation curve is not statistically inconsistent with the observed data. 

Letting p be the underlying actuation probability for some range bin, p (a point 
estimate for p) can be determined by dividing the number of actuations observed in that 
range bin (x) by the total number of observations in that bin (n). We intend to calculate an 
approximate 100(1-a)% confidence region about fp. The confidence region is defined by p, 
and p,,, where: 


- pr. satisfies P(X x)= a/2, where X is a binomial random variable with 
parameters p, and n. ' 


- Py satisfies P(X < x)= «/2, where X is a binomial random variable with 
parameters p,, and n. 
Expressing the binomial distribution explicitly, 


- pf, satisfies > " py (- py" = «/2 


- py satisfies " pu (- py)" = &/2 
i= 0 














If the true binomial actuation probability in a range bin is f,, then the probability of 
receiving the number of actuations actually observed or more is &/2. That is, Pp, 1s taken to 
be the smallest actuation probability which is reasonably consistent with the observed data. 
Sunilarly, if the true binomial actuation probability in a range bin is Py, then the probability 
of receiving the number of actuations observed or fewer is «/2. In this case Py is interpreted 
as the largest p reasonably consistent with the observed data. 

For a= .1 and .2, Figures 4a and 4b show the confidence regions for each range bin 
of the asymmetrical actuation curve previously presented, with an overlay of the symmetric 
actuation curve of Figure 3a. It can be seen that the symmetrical actuation curve is well 
within the boundaries of the confidence regions for these a values. Based on these 
confidence regions and in spite of the asymmetry of Figure 3b, we conclude that it is not 
unreasonable to assume that the actuation curve which generated this data is symmetric. 

2: STANAG 1142 Estimate of Actuation Width and Actuation Probability 

Once an empirical actuation curve has been constructed, the first step delineated in 
STANAG 1142 is to determine the aggregate actuation width (W) by calculating the area 
underneath the empirical actuation curve. The next step in the procedure is to find that 
distance y, which contains the central two thirds of the total area under the actuation curve 
(refer to Figure 3). The actuation probability (B) is then assumed to be the average value of 
the actuation curve over the central two thirds region, and the actuation width (A) is selected 


so that A * B= W. When applied to the actuation data analyzed here, STANAG 1142 
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Figure 4a. Symmetrical Actuation Curve Compared 
with Confidence Region Derived from an 
Asymmetrical Actuation Curve (a=.1). 
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with Confidence Region Derived from an 
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produced estimates for A and B of 1836 yards and .3432 respectively for the symmetric - 
actuation curve of Figure 3a, and 1276.5 yards and .3909 for the asymmetric curve of Figure 


3b. 


D. MAXIMUM LIKELIHOOD APPROACH 

For each actuation curve type, we used maximum likelihood estimates to find those 
parameters resulting in the greatest probability of receiving the actuation results actually 
observed during the exercise. But before proceeding it is necessary to estimate the 
minesweeper’s navigational error. 


1. Standard Deviation and Mean of the Ship’s Navigational Error as 
Determined by Exercise Data 


Two of the parameters needed for the likelihood analysis are an estimate of the 
standard deviation (a) and mean (w) of the ship’s navigational error. For computational 
convenience, these two parameters were estimated by conducting a separate analysis. An 
examination of the exercise’s physical setting suggested a simple approach for estimating 
these parameters. Here, reported range is defined as the lateral distance between the ship 
and a mine as determined through post-exercise reconstruction. Similarly, predicted range 
is that same distance as calculated by assuming the ship followed its intended track exactly. 
Figure 5 is a graphical representation of reported ranges versus the predicted ranges. A 
histogram of the difference between reported range and predicted range was generated, and 
a normal curve was fitted to it, shown in Figure 6. The resulting normal curve yielded 
estimates for the standard deviation and the mean of 159.2 yards and 19.74 yards 


respectively, with a Chi-Square value of 3.69276 which indicates a relatively good fit. 
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Figure 5. Reported Range versus Predicted Range. 
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2. Maximum Likelihood Analysis Background 

An important part of the likelihood analysis was to compute the probability of an 
actuation given a predicted lateral range d. The eventual likelihood function was the product 
_ of these probabilities, one for each ship-mine encounter. The initial assumptions made were: 
normal navigational errors with standard deviation go and mean yw, a symmetrical rectangular 


actuation curve with parameters A and B, and a predicted lateral range of d Letting, 


P(mine actuates | lateral range < A/2, predicted lateral range= d)= B and 


P(lateral range < A/2 | predicted lateral range= d) = fae) = ofA) 
we have 
P(mine actuates | predicted lateral range= d) = 
pe [of2A2H) af e-AR Hy 


Figure 7 is a visual representation of Equation 1. In this figure, the Gaussian curve 
represents the distribution of the ship’s lateral range. The rectangular curve is a generic 
actuation curve for a mine. The probability of actuation is B times the area under the density 
curve within 4/2 of the mine position. 

3. Maximum Likelihood Estimate for Parameter B 

The analytical form of Equation 1 allows us to quickly check if the maximum 
likelihood estimate for B is equal to one. The likelihood function is: 


I 
B's]Ty,, J=0 
r=] 


L(B)= [2] 


I mf 
BT Ly, + [[a-B+y), vet 
i=] j=l 
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where / is the number of positive actuations, J is the number of negative actuations, and y, 
is the probability that the ship comes within A/2 of the mine during mine encounter k. 
Equation 2 is the probability of receiving the actuation results actually obtained from the 
exercise, assuming probabilistic independence for each mine encounter. In order to 
determine the maximizing B, the log of the likelihood function was differentiated with 


respect to B and set equal to zero. 


el. 2. 5) 
1 ‘1 


[3] 
dB B j= -B *y 


which gives the following when evaluated at B=1: 


AlogiL(B)ly _ op _ 
dB ea >| 





ed [4] 


Due to the concavity of log [Z(B)] for 0< B <1, Equation 4 implies that if the number of 


J 
actuations I > > (+), then the slope of log [Z(B)] is positive at B=1 and the maximizing 
j=l m4 j 





B (in the interval 0< B <1) must be one. IfJ< > ( = ), then the maximizing B is found 
by numerically solving for the root of Equation 3 in the oawal O< B <1. The importance 
of this procedure is that it allows us to solve for the best B once the other model parameters 
have been specified. This reduces the amount of searching we have to do to find the model 
parameters maximizing the likelihood function. 

4. Results of the Analysis 


Using the theoretical background previously discussed, three MATLAB models were 


developed. The models were used to evaluate the exercise data and to compute those 
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Figure 7. Probability of Actuation. 


parameters corresponding to the greatest likelihood of receiving the actuation results actually 
obtained, had the data come from the particular actuation curve being modeled. Ali three 
models were very similar, with the major difference being in the expressions used to 
compute the conditional probability of a mine actuation (see Equation 1). 

a. Symmetrical Rectangular Curve 

The first and simplest case considered was a symmetrical rectangular 
actuation curve. In this model, Equation 1 remained unchanged from the form presented 
during the discussion of the theory. For this case the parameters of interest were actuation 
width (4), and actuation probability (B). The MATLAB code for this model is included as 


Appendix B, its run time was 75 seconds on a 486/33MHz PC. 
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b. Asymmetrical Rectangular Curve 
The second model examined the fit of the data to an asymmetrical rectangular 


actuation range curve. Here Equation 1 became: 


P(actuation | predicted lateral range= d) = 


pe [SAH ofA] 





In this model the parameters of interest were actuation width ( as defined by A, and A, in 
Figure 2) and actuation probability (B). The MATLAB code for this model is included as 
Appendix C. The run time for this model was 11 minutes . 

c Washburn’s Curve 

The final case examined the use of a more complicated, yet more realistic, 
type of actuation curve to analyze the exercise data. Using Washburn’s actuation curve the 


conditional probability of a mine actuation was computed as: 


P(actuation | predicted lateral range = d) = 


I-A ew (2A, 


co 


[le +e 


- A o+/In 





where x is the lateral range; and A, B, and c are parameters. The MATLAB code for this 
model is included as Appendix D. The run time for this model was 39 minutes. 

As previously discussed, the exercise data is assumed to have been generated from 
the symmetrical actuation data shown in Figure 3a. However, for the purpose of illustrating 


the use of the methodology for asymmetric data, the three models were also run under the 
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assumption that the data had been generated from the asymmetrical actuation data of Figure 


3b. The results for both data sets are summarized below in Tables | and 2: 


Rectangular Symmetrical Rectangular Asymmetrical Washburn Symmetrical 
Actuation Curve Actuation Curve Actuation Curve 













Fe 


Table 1. Summary of Results Using Symmetric Data. 


Table 1 shows the symmetric data results. As expected, the rectangular asymmetric 
and the rectangular symmetric fits were essentially identical. What was more surprising was 
that the best Washburn fit was also rectangular. This is shown by the large value of the best 
fit C. Experimentation with the data revealed that more actuations with small lateral ranges 


increases the central tendency of the data and depresses the best fit C. 


Rectangular Symmetrical Rectangular Asymmetrical Washburn Symmetrical 
Actuation Curve Actuation Curve Actuation Curve 
A (Yds) 2325 1175 (A,= 100, A,= -1075) 2350 
4286 3376 


Ce 
Likelihood | 2.24 x 10 5.38 x 10” 2.24x 10 


Table 2. Summary of Results Using Asymmetric Data. 
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When the three idealized actuation curves were fit to the asymmetric data, the 
rectangular asymmetrical actuation curve provided the best fit, as was expected. Once again, 


the best fit Washburn curve closely approximated the rectangular symmetrical fit. 





ae 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


This thesis develops a methodology for the evaluation of mine actuation data. It 
makes use of maximum likelihood estimates (MLE) theory to determine the desired 
parameters, for various mine actuation curves, corresponding to the largest likelihood of 
obtaining the observed data. These parameters are important in mine warfare because they 
are inputs to larger models used to plan and evaluate mine clearance operations. The main 
advantages of this methodology are its applicability to both symmetric and asymmetric data, 
and its conceptual simplicity. Its main drawback is the large computational effort required 
to maximize the log likelihood function. Separate analysis methods are also presented for 
testing the symmetry/asymmetry of exercise data, and for estimating the mean and standard 
deviation of a ship’s navigational error. 

The methodology presented here has a direct application to mine warfare and could 
be of considerable use to minefield planners and analysts. As such, further research should 
be conducted to expand it and exploit fe potential benefits. Possible areas of future work 
include: 

@ Applying the methodology to larger data sets (both symmetric and 

asymmetric). 

® Incorporating the methodology into existing tactical decision-aids. 

@ Developing a graphical user interface to facilitate the use of this methodology 

in the Fleet. 


@ Exploring the use of logistic regression to fit the empirical actuation data. 
[Ref. 11] 
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APPENDIX A. EXERCISE MINE ACTUATION DATA 


Number 


9 (North) OR150 XL8 0 OCO00 XL540 XL510 
10 (South) | 0L1340 OLLO90 XL1060 OL495 AL440 
















O- Non-actuation 

X- Actuation 

R- Right 

L- Left 

C- Centerline 

Numbers- Lateral distance in yards 


Note: In the table above, the result of the mine encounter 
during track 1 for mine 7 was: no actuation (0), the encounter 
occurred on the right in relation to the ship’s track (R), 

and the lateral distance between the mine and the ship was 690 


yards (690). 
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APPENDIX B. MATLAB CODE FOR SYMMETRICAL RECTANGULAR 
CURVE USING SYMMETRIC DATA 


clear 
cle 


Sigma= 159.2; 

mu= 19.74; 

d= [690,410,410,190,205,1540,1160,1170,675,680,550,275,... 
235430 ,3550 7410059507990) 3354297,280,2004 60 ;-3005-4.8% 
18350;900,640,°700,160,0;,,90,210,340,900,1030; 780) 4.%- 
54:0;,:560.,:25,4°75; 150; 80 -0..540,510,1340;,1090,1060,475....2« 
440]; 


d neg= [690,410,410,190,1540,1160,1170,675,550,275,235,... 
35071200, 335; 160,180,800, 1850; 900,640,700, 1605 4.- 
90 ;340;,,900,1:030;,780,540, 560, 75,150 ;0,1340,1090:. 4 «4 


475]; 
Global: “yy -neg; 
Pval= 0: | 
min Pval= 1000; | | . | 
A= 0; 
B= 0 
i. aCl= ~L5% 





FOr T= 1:2 200+ 
a= O+ (1-1)*25; 
y_neg= snormcum((d_neg+(a/2)-mu)/sigma)-... 
snormcum( (d_neg- (a/2) -mu) /sigma) ; 
n= sum(y_ neg./(1-y_ neg) ) ; 
LE nNoact Sn 


biti) = 1 
else 
D(a )= Pzero(" ft", ..35) = 
end 
Prob=b (1) * [snormcum( (d+ (a/2)-mu) /sigma) -... 
snormcum( (d- (a/2) -mu) /sigma) ] ; 
Pd= (det_out.* Prob) + [(1-det out) .*(1-Prob)]; 
Pval(1)= sum(abs(log(P 3 


1f min Pval > Pval(i), 


min Pval = Pval (i); 
A= a; 
B= b(1); 

end 


end 





forintf(’The likelihood is= %*g.\n\n’, exp(-min_Pval) ) 
forint’ (A= sg. \n\n" ,..A) 
forinte( B= -2¢q.\i\ ni’. -B) 
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APPENDIX C. MATLAB CODE FOR ASYMMETRICAL RECTANGULAR 
CURVE USING ASYMMETRIC DATA 


clear 
cle 


sigma= 159.2; 

mu= 19.74; 

d= [690,410,410,-190,-205,-1540,-1160,-1170,~-675,-680,... 
SIU 61 Op 29 O¢ MSL Oy 350, A 11005950 7-990 7H 335) S20 ed 
180,-100,-180, -800,-1850, -900, -640,-700,-160,0,90,... 
=21.0;, +340,~2900,-1030,-780,-540, =560 25,75; 150;)-80,07 ..% 
-540,-510,-1340,-1090, -1060, -475, -440] ; 


d neg= [690,410,410,-190,-1540,-1160,-1170,-675,550,275,... 
23543 00721100; =3357180;>-160,+800;.-1650;-900)... 
#640 ,--700, 526 0,.90,2340,-900,.=1030,-780;,-540,-560;.42% 
75,150,0,-1340,-1090,-475] ; 


global. ‘y neg; 


Pyval= 0; 

min. Pval= 1000; 
A= Q; 

B=: 02 


n act= 15; 


FOr a= 1:44; 

al= 04 (ie1)*25; 

LOL FS 2243; 
a2= 1050- (3-1) *25; 
y_neg= snormcum((d_neg+al)-mu/sigma) -... 

snormcum((d neg-a2) -mu/sigma) ; 

n= sum(y_neg./(1-y neg)) ; 
Lt sy act Sn 


Died) 1S 
else 
, tiga) = t7erel* o42.5)3 
end 


Prob=b (1,3) * [snormcum( (d+a1)/sigma)-... 
snormcum ((d-a2) /sigma) ] ; 

Pd= (det_out.* Prob) + [(1-det out) .*(1-Prob)]; 

Pval(i,j)= sum(abs(log(Pd))); 

1£ min Pval > Pval(i,j), 


min Pval = Pval(i,j); 
Al= gai; 
A2= a2; 


35 


end 
end 
end 


fprintf 


EPrInGE 
EprInce 


COTA 
Poriner( - 7. 
('A2 
(‘B= 


e likelihood is= %g.\n\n’, exp(-min_Pval) ) 
2 2G.\nn > Aly 

= %g.\n\n’, A2) 

g.\n\n', B) 
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APPENDIX D. MATLAB CODE FOR WASHBURN SYMMETRICAL 
CURVE USING SYMMETRIC DATA 


clear 
cle 


Sigma= 159.2; 
mu= 19.74; 


d= [690,410,410,190,205,1540,1160,1170,675,680,550,275,... 
23973107350 ;,1100,950,5.9907335,297; 1805.100,180-,800. oa: 
1850,900,640,700,160,0,90,210,340,900,1030,780,... 
240; 560925: 15g 150400707 540,,5 1.0, 1340,1090,1060,475).5<. 
440]; 


d neg= [690,410,410,190,1540,1160,1170,675,550,275,235,... 
350,1100 ;335,180,180,800;1850, 900,640; 700, 160,22. 
90,340,900,1030, 780,540,560, 75,150,0,1340,1090,... 
475]; | 


x Neg= -1/91.4375106..5721900; 
xe SLI0O S10 F900 


global y_neg; 


Pval= 0; 

min Pval= 1000; 
A= Q; 

B=. 0; 

C=O: 


Heer 157 


for 1] 4230; 
EPrianer ("= eo .\nn’ ¢, 4) 
a= 500+ (1i-1)*25; 
C= 03 
LOL y= 1.310% 
C= c+ 310+ 
Lor Ha Ls35 
S= d_neg(h) -x_neg; 
a_neg= exp(-abs((x neg/a)).*c); 
b neg= exp(-.5*((s-mu) /sigma) .*2); 
Z neg= a_neg.*b neg; 
y_neg(h)= 1/(sigma*sqrt (2*pi))*simrule(z neg,65.71) ; 
end 
n= sum(y_neg./(1-y_neg)); 
Ee De ack. > 
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b= iis 
else 

b= £zero(* tf"... 5)4 
end 
for k= 1250 

r= d(k) -x; 


z= exp(-abs((x/a)).*c) .*exp(-.5*((r-mu) /sigma) .*2) ; 
Prob(k)= (b/(sigma*sqrt (2*pi)))*simrule(z,10); 

end 

Pd= (det _out.* Prob) + [(1-det_out).*(1-Prob)]; 

Pval= sum(abs(log(Pd))),; 

i1£ min Pval > Pval 


min Pval = Pval; 

A= 2%*a; 

C= a 

B= jb; 

end 
end 

end 
fprintf(’The likelihood is= *g. \n\n’, exp(-min_ Pval) ) 
Porineel (A= sq. \n in’, A) 
Forinur(’ C= 4g. nm’; “C) 
fprintt (’ B= sg.\n\n', B) 
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